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DISCHARGE  OF  FIRE  SUPPRESSION  AGENTS  FROM  A  PRESSURISED  VESSEL; 
A  MATHEMATICAL  MODEL  AND  ITS  APPLICATION  TO  EXPERIMENTAL  DESIGN 
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Fire  Modeling  Group 
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ABSTRACT 

A  mathematical  model  and  associated  computer  program  is  developed  to  simulate  the  discharge  of 
fire  extinguishment  agents  from  Nj-pressurized  vessels,  ihe  model  is  expected  to  have  three 
applications.  First,  to  establish  an  experimental  design  and  procedure  which  closely  simulates 
discharge  of  a  field-deployed  vessel;  second,  to  evaluate  the  discharge  characteristics  of  a  wide  range 
of  altemative-agent/pressure-vessel  configurations,  thereby  extending  the  slow  and  relatively  costly 
experimental  method  of  making  such  evaluations;  and  finally,  to  predict  vessel  exit  flow  conditions 
to  be  used  to  solve  the  problem  of  agent  dispersal  outside  of  the  discharge  vessel.  The  model  is  used 
in  example  calculations  which  address  the  first  of  these  applications.  . 

The  field-deployed  system,  which  forms  the  basis  of  the  example  calculations,  involves  a  half-liter 
cylindrical  discharge  vessel  with  a  circular  discharge  nozzle/orifice  of  diameter  0.019m.  The  vessel 
is  half-filled  with  liquid  Freon  22  and  is  pressurized  with  N2  to  41.37xl0sPa  (600psi).  Vessel 
discharge  is  initiated  by  actuation  of  an  explosive  cap  over  the  nozzle/orifice. 

The  simulating  experimental  configuration  involves  a  modified  field-deployed  system.  A  diaphragm 
with  nominal  41.37xlOsPa  (600psi)  rupture  pressure  [actual  values  between  37.92xl05Pa  (55Gpsi)  and 
44.82xl05Pa  (650psi)]  replaces  the  explosive  cap.  The  system  is  equipped  with  a  high-pressure  N2 
holding  tank  connected  to  the  discharge  vessel  via  an  orifice.  An  experimental  run  begins  with  the 
onset  of  through-orifice  N2  flow  from  the  holding  tank.  The  vessel  is  pressurized  to  the  point  of 
diaphragm  rupture  and  this  is  immediately  followed  by  vessel  discharge. 

The  model  is  used  to  simulate  discharge  of  the  field-deployed  system  and  pressurization/discharge  of 
the  experimental  system.  Simulations  of  the  experimental  system  involve  holding  tank  volumes  of 
23x1 0'3m3  or  2.5xl0‘sm3;  orifice  diameters  of 0.005m,  0.001m,  or  0.0005m;  and  initial  vessel  pressures 
of  938x1 O^Pa  (136psi)  (the  saturation  pressure  of  Freon  22  at  294K)  and  34.47xl05Pa  (500psi). 

From  the  calculations  it  was  determined  that  the  2.5xl0‘3m3  holding  tank  with  the  0.0005m  orifice 
could  be  used  to  simulate  accurately  the  discharge  of  the  field-deployed  system  and  that  it  is 
reasonable  to  expect  that  this  experimental  design  would  give  good  simulations  even  when  extended 
to  a  range  of  parameters  and  agent  materials  well  beyond  the  scope  of  the  present  calculations. 
Calculations  also  indicated  that  use  of  the  23xl0‘5m3  holding  tank  and/or  the  0.005m  orifice  would 
not  be  consistent  with  an  acceptable  experimental  design. 

Keywords:  agents  aircraft  fire  safety,  discharge,  fire  extinguishment,  fire  safety,  baton,  halon 

alternatives 


THE  PROBLEM  AND  THE  OBJECTIVE 


This  work  formulates  a  mathematical  model  to  simulate  the  discharge  of  Halon  and  Halon-alternative 
fire  extinguishment  agents  from  Nv pressurized  vessels.  The  objective  is  to  develop  a  mathematical 
modci  which  simulates  agent-discharge  experiments  now  under  way  at  the  Building  and  Fire  Research 
Laboratory  (BFRL)  of  the  Nat  onal  Institute  of  Standards  and  Technology  (NIST).  The  experiments 
are  part  of  a  program  to  support  advances  in  fire  safety  in  aircraft  of  the  US  Air  Force. 

The  model  is  expected  to  have  three  applications.  First,  it  will  be  used  to  establish  an  experimental 
design  and  procedure  which  closely  simulates  discharge  of  field-deployed  vessels  while  allowing  for 
acquisition  of  data,  including  high  speed  photography,  to  characterize  adequately  the  discharge 
process.  Second,  the  mode!  will  be  used  to  evaluate  the  discharge  characteristics  of  a  wide  range  of 
alternative-agent,  pressure-vessel  configurations,  thereby  extending  the  slew  and  relatively  costly 
experimental  method  of  making  such  evaluations.  Finally,  it  will  be  used  to  determine  the  discharge 
vessel  exit-flow  conditions  Lor  use  in  the  simulation  of  agent  dispersal  outside  of  the  vessel.  .After 
presenting  the  mathematical  model,  this  w*ork  will  include  example  calculations  which  address  the  first 
of  these  applications. 

The  analysis  is  based  on  the  experimental  arrangement  depicted  in  Figure  1. 

THE  EXPERIMENT  AND  THE  MODEL  ASSUMPTIONS 
The  Experimental  Arrangement 

Refer  to  Figure  1.  This  represents  the  arrangement  associated  with  discharge  experiments  being 
considered  at  NIST  The  Figure  1  arrangement  can  also  be  used  to  describe  the  phenomena  of 
discharge  from  vessels  under  field  conditions. 

The  experimental  arrangement  involves  a  right  cylinder  discharge  vessel  at  pressure  PqV,  which  would 
contain  N2-pressurized  test  agent,  and  a  holding  tank  filled  with  mass  of  N2  at  pressure  Pj^j. 

The  discharge  vessel  is  of  height  ZDV  and  cross-sectional  area  ADV.  The  volume  of  the  holding  tank 
is  Vm. 

The  holding  tank  is  connected  via  an  orifice  of  area  Aq  to  the  discharge  vesseL  It  is  assumed  that 
the  flow  path  through  the  orifice  can  be  opened  or  closed  with  a  relatively  fast-acting  solenoid  valve. 
If  the  path  is  open  it  is  assumed  that  Pj^  >  PDV,  i.c.,  the  flow  is  always  from  the  holding  tank  to  the 
discharge  vessel.  In  the  case  of  a  field-deployed  system  there  is  no  holding  tank  and  Aq  is  zero. 

At  the  bottom  of  the  vessel  is  a  shor;-nozzle/orifice-type  opening  of  area  AN.  The  agent  will  be 
discharged  through  this  to  the  outside  ambient  environment  which  is  at  pressure  P^mb  <  PDV.  The 
discharge  flow  path  is  originally  closed  off  by  a  cap  or  diaphragm.  When  the  cap  is  removed,  the  test 
agent  liquid  is  driven  out  of  the  discharge  vessel  by  virtue  of  the  cross-nozzle/orifice  pressure 
difference. 
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The  Procedure  Prior  to  An  Experimental  Run  or  a  Field-Deployed  Discharge 

The  orifice  and  nozzle/orifice  flow  paths  are  closed  and  the  discharge  vessel  is  evacuated.  The  vessel 
is  then  filled  completely  with  a  known  mass  of  test  agent.  There  will  be  a  volume  of  liquid  agent 
below  and  a  volume  of  gaseous  agent  above  and  the  pressure  will  be  PSAT,  the  saturation  pressure 
at  the  agent  temperature.  The  vessel  is  then  pressurized  with  N2.  This  flows  into  the  vessel  from 
the  holding  tank  or  from  some  other  relatively  high-pressure  N2  source.  In  general,  the  upper  gas 
volume  is  now  a  two-component  mixture  of  N2  and  test  agent  gas.  Although  some  N2  may  be 
dissolved  in  the  liquid  volume,  it  is  assumed  that  the  amount  is  always  so  small  that  throughout 
subsequent  sequential  (additional-)pressurization  and  discharge  processes  to  be  studied  here  the 
properties  of  the  liquid  in  the  discharge  vessel  are  well  approximated  by  the  liquid  properties  of  the 
pure  test  agent. 

Now  consider  the  system  at  the  time,  t  =  0,  when  an  experimer  al  run  or  a  field-deployed  discharge 
is  initiated.  The  liquid/gas  interface  is  a  distance  Zx  above  thv.  bottom  of  the  vessel.  The  mass  of 
liquid  in  the  vessel  is  M^val  subscript  refers  to  that  portion  of  Discharge  Vessel  test  Agent  in 
the  Liquid  state)  and  throughout  the  experimental  run  or  Geld-deployed  discharge  the  temperature 
and  density  of  the  liquid  that  remains  in  the  vessel  will  be  assumed  to  remain  at  the  now-existing 
respective  values  Tal  and  =  PsaT<Tal)-  As  indicated,  is  estimated  to  be  the  density  of 
saturated  test  agent  liquid  at  the  saturation  temperature  Tal-  The  gas  mixture  is  at  an  assumed 
uniform  initial  temperature  TDV,l  and  the  initial  masses  of  gasec  xs  test  agent  and  N2  in  the  discharge 
vessel  are  MDVj\g  anc*  ^dvnz.1’  respectively.  Note  that  the  re-pressurization  process,  prior  to  t 
=  0.  may  have  been  so  rapid  that  the  liquid  and  gas  are  not  in  thermal  equilibrium  with  each  other. 
For  this  reason,  TDV1  is  not  necet:arilv  identical  to  Tal-  The  discharge  vessel  and  the  holding  tank 
are  at  the  uniform  pressures  PDV1  and  respectively. 


During  an  Experimental  Run  or  Field-Deployed  Discharge 

During  discharge  of  a  Figure  1-type  system,  the  flow  path  through  the  orifice  is  either  open  or  closed 
at  t  =  0+.  When  the  orifice  is  closed  the  discharge  will  simulate  discharge  of  a  field-deployed  system 
which  does  not  involve  a  holding  tank.  When  the  orifice  is  open  the  discharge  will  simulate  discharge 
of  the  experimental  system  where  additional  N2  pressurization  during  the  discharge  is  provided  by  the 
holding  tank.  The  flow  path  through  the  nozzle/orifice  at  the  exit  of  the  discharge  vessel  is  either 
opened  at  t  =  0+  (as  in  the  field-deployed  system,  which  uses  an  explosive  cap  to  initiate  discharge 
from  an  equilibrium  state  at  the  specified  PDV1)  or  e^e  opens  if  and  when  PDV  rises  to  some 
specified  diaphragm  rupture/burst  pressure,  PguRST  <  ^hti  (“  'n  the  experimental  system).  The 
time  when  this  latter  flow  path  is  opened  is  designated  as  tBURST- 


Further  Assumptions 

For  times  and  conditions  of  interest  it  is  assumed  that  as  long  as  liquid  remains  in  the  vessel  PDV  > 
Psat(Tal)  and  there  is  no  possibility  of  flashing  (i.e.,  spontaneous  change  from  a  thermodynamically 
unstable  liquid  state  to  an  equilibrium,  two-phase,  liquid/gas  state)  in  the  liquid  volume.  Future  work 
will  address  the  problem  of  removing  this  latter  constraint  from  the  analysis.  It  is  assumed  that  any 
N,  dissolved  in  the  liquid  that  may  come  out  of  solution  during  the  discharge  process  is  negligible. 
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The  purpose  *.'f  this  work  is  to  mathematically  model  the  state  of  the  system  at  any  t  >  0  up  to  tD, 
where  tD  is  the  smaller  of  1)  the  time  when  the  discharge  vessel  is  emptied  of  test-agent  liquid,  and 
2)  the  time  when  PDV  is  reduced  to  P$at(^al)'  ^°r  experimental  systems  of  interest  here  tgURST 
is  expected  to  be  of  the  order  of  10s  and  the  time  interval  of  the  discharge  process,  (tD  -  tBURST), 
of  the  order  of  l(T~s.  Prior  to  tD  there  is  no  gas  mixture  discharge  from  the  vessel.  The  mass  of 
gaseous  agent  in  the  vessel  therefore  remains  constant  at  its  initial  value  MDV^G. 

In  view  of  the  relatively  short  times  of  interest  it  is  assumed  that  during  an  experimental  run  or  field- 
deployed  discharge  there  is  no  heat  or  mass  transfer  across  the  liquid/gas  interface.  It  is  also  assumed 
that  there  is  no  heat  or  mass  transfer  across  the  interface  that  contains  the  nitrogen  and  gaseous  test 
agent  in  the  combined  holding-tank/discharge -vessel  system,  e.g.,  there  is  negligible  heat  transfer  to 
the  walls  of  the  vessel  and  holding  tank.  This  is  the  basis  of  the  assumption  that  both  MDVtAG  and 
the  total  mass  of  N2  in  the  combined  system,  designated  as  MN2,  are  constant  for  all  t,  0  <  t  <  tD. 
MN2  is  determined  from  the  initial  conditions. 

It  is  assumed  that  the  gaseous  test  agent  and  the  N2  can  be  modeled  as  perfect  gases  with  constant 
specific  heats  and  gas  constants,  CV^G  Cp^G,  RAG  and  CVN2  CPN2,  RN2.  The  values  of  the  specific 
heats,  which  depend  on  temperature,  are  taken  to  be  those  which  correspond  to  tdv,i- 

THE  MODEL  EQUATIONS 

Together  with  specified  parameters,  the  unknown  time-dependent  variables  of  the  problem  that 
define  the  state  of  system  at  an  arbitrary  time  are:  Pj-rr*  Ppv»  ^HTN2*  ^DVN2*  ^dv,al*  ^HT*  *Pdv* 
Z.  Known/specified  parameters  would  include:  the  initial  values  of  these  variables,  Phti>  Pdvi* 
Mht,N2,i*  ^DVN2,i’  ^HXi’  Tdv,i*  anc*  respectively;  the  geometric  parameters  of  the 

system,  Aq,  An,  Adv,  Zdv,  V^;  material  properties  of  the  test  agent  and  of  N2;  the  constant  value 
^DVag;  an<*  nozzle/orifice  discharge  coefficients,  to  be  introduced  below. 

The  equations  governing  the  variables  are: 


PHT  -  (MHXN2^VHt)RN2THT  0) 

?hV  =  (Mdv^gRag  +  MDV^2RN2)TDV/[(ZDV  *  Z)Adv]  (2) 

Conservation  of  Mass  of  N-, 

Mdv,N2  +  Mjj-£N2  =  constant  =  MN2  =  MDVrN? 1  +  (3) 
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Relation  of  Z  to  Mpy^u. 

^DVjAl/^AL  =  AdV^  W 

Reversible  Adiabatic  Expansion  for  N2  in  Holding  Tank 

^m/(^HXN2)y  =  constant  =  PHT;i/(MHT^2,i)y  @) 

where  y  =  v,N2  =  ^N2^yN2  +  1- 

First  Law  of  Thermodynamics  for  Entire  Gas  System 

d[(MDVAG^yAG  +  MDyN2CV,N2)TDV  +  MHXN2CyN2THT]^dt 

=  rate  of  work  done  on  the  gas  system 
=  PDVADVdZ/dt  (6) 


Flow  Across  the  Orifice 


-  cD>0AoPHT{[W(RN2Tm)](2/(y  +  i)](y  +  1V(y  •  x))xa 


d^HXN2^dt  ~ 


ifPm^DV  *  [(y  +  l)/2]ty/(y'1)1  (Le.,  choked  flow) 

(7) 


-  Cd  oAqPj^Ppv/Pht)1^^/!!  -  (P DV^P HT)^y  '  X)^(y  *  1)^2^^]^ 

if  Pht/Pdv  <  [(y  +  l)/2]^,(ym1^  (i.e.,  un-choked  flow) 


where  CD  0  is  the  (compressible)  flow  coefficient  for  the  orifice  [1]. 


Flow  across  the  Nozzle/Orifice 


dMoYAi/dt  - 


*  CD>nAN[^al(PdV  *  PAMb)]1/2  **”  PAMB  -  PSP 

*  Cd,nAn[^al(pdv  *  Psp))1/2  PAMB  <  PSP 


(8) 


Where  Cp  ^  is  the  flow  coefficient  for  the  nozzle/orifice  and  P<jP  is  defined  below.  Assuming  that 
An/Adv  <  <  1  in  discharge  vessels  of  interest  here,  in  Eq.  (8)  the  kinetic  energy  of  the  liquid 
upstream  of  the  exit  nozzle  is  neglected  compared  to  the  kinetic  energy  at,  and  immediately 
downstream  of  the  nozzle. 


DEFINITION  OF  Psp  AND  COMMENTS  ON  EQ.  (8) 

In  Eq.  (8)  it  is  assumed  that  the  liquid  in  the  discharge  vessel,  at  state  [XP]  =  [T^,  PDV  > 
psatOal)]’  flows  into  and  through  the  exit  nozzle/orifice  while  moving  along  state  paths  of  constant 
entropy.  When  the  pressure  of  the  liquid  drops  below  its  saturation  pressure,  the  liquid  is  assumed 
to  initiate  its  movement  into  the  "vapor  dome"  as  a  metastable  super-heated  liquid. 

For  the  generic  pure  material,  the  relevant  P-V  diagram  which  depicts  the  metastable  states  is 
sketched  in  Figure  2.  This  includes  the  region  of  metastable  super-heated  liquid  states  and  the 
corresponding  region  of  metastable  subcooled  vapor  states.  For  a  pure  material  the  metastable  liquid 
or  vapor  state  can  only  be  maintained  where  (3P/3V)  |T  <  0  [2],  Le.,  only  outside  the  locus  of  points, 
referred  to  as  the  spinoidal  curve,  where  (3P/3V)  |T  =  0.  The  spinoidal  curve,  sketched  in  Figure  2, 
passes  through  the  critical  state  and  has  the  superheated  liquid  region  to  its  left  and  the  subcooled 
vapor  region  to  its  right. 

At  any  time  during  the  discharge,  Psp  is  defined  as  the  particular  pressure  along  the  liquid-leg  of  the 
spinoidal  curve  of  the  test  agent  associated  with  the  intersection  of  the  spinoidal  curve  and  the  above- 
mentioned,  instantaneous,  constant-entropy  paths.  Since  states  along  the  spinoidal  curve  are  states 
of  unstable  equilibrium  where  spontaneous  nucleation  will  occur  [2,  3],  Le.,  the  constant-entropy 
metastable-state  path  of  the  liquid  cannot  be  sustained  at  the  state  defined  by  the  intersection  point, 
it  is  conjectured  that  if  PSP  >  P violent  flashing  of  the  liquid  jet  occurs  where  and  when  its 
pressure  reaches  Psp.  This  would  be  upstream  of  a  physically-unachievable,  Le^  unstable,  vena 
contractu  at  pressure  P**®. 

ff  pamb  **  achieved  prior  to  the  P§p  state  then  the  top  portion  of  Eq.  (8)  is  used.  For  this 
circumstance  the  constant-entropy  state  of  the  possibly  metastable  superheated  liquid  at  P^^  is  a 
physically-admissible  endpoint  to  the  liquid-jet  development  process.  The  situation  is  consistent  with 
the  idea  that:  1)  the  endpoint  state  coincides  with  the  vena  contracts  of  the  jet;  2)  the  jet  is 
convected  for  some  distance  downstream  of  this  endpoint  without  substantial  changes  and  at  the  vena 
contracta  diameter;  and  3)  the  jet  eventually  breaks  apart  due  to  fluid-dynamic  and/or  thermodynamic 
(flashing)  instabilities.  Note  that  the  upper  Eq.  (8)-description  of  the  nozzle  flow  rate  of  a 


superheated  liquid  as  an  incompressible,  non-flashing  fluid,  is  consistent  with  results  reported  in  the 
literature,  e.g.,  [4],  for  flow  through  sharp-edge  orifices  and  short  nozzles. 

If  PSp  is  achieved  prior  to  P^mb  then  flashing  of  the  metastable  liquid  will  be  initiated  immediately 
downstream  of  the  position  in  space  where  the  intercept  occurs.  In  the  latter  case  the  incompressible 
flow  calculation  methodology  is  still  applicable  at  and  upstream  of  the  spinoidal-curve  intercept  point. 
This  is  indicated  by  the  use  of  the  bottom  portion  of  Eq.  (8). 

The  above  discussion  is  illustrated  with  C02.  For  stable  and  metastable  liquid  C02,  a  sketch  from 
[3]  of  constant-entropy  paths  in  a  P-T  diagram  is  presented  in  Figure  3.  In  preparing  this  figure,  the 
Peng-Robinson  equation  of  state  was  used  in  [3]  to  describe  the  metastable  liquid  C02. 

First  assume  that  liquid  C02  is  in  a  Figure  1-type  discharge  vessel  at  PDV  =  TOxK^Pa  and  T^  = 
288K.  In  Figure  3  this  initial  state  is  seen  to  lie  on  the  s  =  s2  constant-entropy  line,  which  has  no 
positive  P  intercept  with  the  spinoidal  curve.  As  the  liquid  CO,  flows  toward  and  then  out  of  the 
vessel’s  nozzle/orifice  its  state  is  assumed  to  move  downward  along  the  s2  curve  of  Figure  3.  As  can 
be  seen,  the  P  =  PATM  =  lxlC^Pa  intercept  on  this  curve,  corresponding  to  T  *  277K,  represents 
an  achievable  metastable  liquid  state  for  the  material.  This  is  the  metastable  liquid  state  that  is 
predicted  by  the  present  model.  It  would  be  expected  at,  and  for  some  distance  downstream  of,  the 
jet’s  vena  contracta. 

Now  assume  that  liquid  C02  in  a  discharge  vessel  is  at  P  =  PDV  =  tjQxK^Pa  and  T  =  T^  =  291 K. 
In  Figure  3  this  initial  state  is  seen  to  lie  on  the  s  =  s3  constant-entropy  line,  which  intersects  the 
spincidal  curve  at  approximately  Psp  =  lOxl^Pa  >  >  P^g.  In  this  case,  as  the  liquid  C02  flows 
through  and  out  of  the  vessel  nozzle/orifice  the  model  predicts  that  it  will  flash  explosively  when  P 
=  Psp.  This  will  occur  upstream  of  the  position  at  which  a  fluid  jet  vena  contracta  would  otherwise 
occur,  well  within  a  distance  of  one  nozzle/orifice  diameter  downstream  of  the  vessel  exit 

This  nozzle/orifice  flow  model  is  consistent  with  a  discharge  process  involving  a  smoothly  time-varying 
PDV,  a  constant  T^,  and  an  abrupt  change  in  nozzle/orifice  exit  flow  from  a  simple  incompressible- 
fluid-jet  type  flow  to  an  explosively  flashing  two-phase  flow.  This  kind  of  discharge  phenomenon  was 
clearly  observed  during  recent  Freon  22  discharge  experiments  at  NIST  [5].  Note  that  the  idea  of 
using  the  condition  of  the  existence  of  Psp  >  PAMB  as  a  criterion  for  violent  flashing  of  liquids 
flowing  through  sharp-edged  orifices  or  very  short  nozzle-like  openings  does  not  seem  to  have  been 
proposed  previously.  Nor  has  it  been  validated  quantitatively. 

REDUCED  DIMENSIONLESS  EQUATION  SET  -  THE  INITIAL  VALUE  PROBLEM 

Eqs.  (1)  -  (8)  are  made  dimensionless  by  introducing  a  dimensionless  time,  r,  dimensionless  mass  of 
N2  and  liquid  agent  in  the  discharge  vessel,  x2  and  x3,  respectively,  and  a  combined  variable,  Xj. 

*1  =  [(1  +  ^3X2)/(1  +  dv/P DVlX^^i  *  *3)  +  ^2*4(1  *  x2? 

x2  ~  MDVN2/MK2!  x3  =  MDV,Al/MDV,AL4  (9) 

M 


T  ~  tCDJtANpPDV,l/(MDyAL,lADvZi)] 


where 


^1  “  (^K2^DVAg)(^-N^Ag)!  ^2  -  0>HXl^>DV,iy(MHXN2,l^N2)>’» 
^3  =  h(YAO  -  my  - 1);  A4  =  [(yAG  - 1  )/(y  -  1)](Vht^(ZiAdv)]; 


(10) 


The  initial  values  of  the  x’s  are  functions  of  the  various  parameters  of  the  particular  problem  of 
interest.  The  dimensionless  initial  value  problem  for  an  arbitrary  Figure- 1  configuration  is 


dx/dT  =  o$  x^O)  specified,  i  =  1,  2,  or  3 


(11) 


where 


ai=(7i(xi»  x2’  xi'  parameters  of  the  problem),  i  =  1,  2,  and  3  (12) 

The  reader  is  referred  to  the  Appendix  for  definitions  of  oi  [the  right  hand  sides  of  the  equations  of 
(A-5)]  and  for  explicit  equations  to  extract  values  of  time-dependent  dimensional  variables  from  a 
solution  of  Eqs.  (11)  and  (12). 


N2  JET-DRIVEN  MIXING  IN  THE  DISCHARGE  VESSEL  -  ESTIMATING  THE  RESULTING 
DISTURBANCE  OF  THE  UQUID/GAS  INTERFACE 

The  above  model  assumes  that  throughout  the  pressurization  and  discharge  process  the  agent  gas  and 
the  N2  in  the  discharge  vessel  are  fully  mixed  and  in  a  state  of  thermodynamic  equilibrium.  It  is  also 
assumed  that  there  is  no  significant  heat  or  mass  transfer  interactions  at  the  gas/liquid  interface, } 
it  is  assumed  that  the  interface  is  relatively  quiescent 

In  cases  where  there  is  no  N2  flow  from  the  holding  tank  (e.g.,  the  orifice  is  closed)  the  initial  fully- 
mixed  state  of  the  gases  will  persist  throughout  de-pressurization,  the  gas  volume  will  be  relatively 
quiescent,  and  the  assumption  of  negligible  interface  interactions  is  expected  to  bold.  This  is  the  case 
for  actual  field-deployed  systems. 

When  there  is  N2  flow  from  the  holding  tank,  it  is  the  N2  jet  from  the  orifice  that  will  drive  gas 
mixing  in  the  discharge  vessel.  Refer  to  Figure  4.  If  jet  velocities  are  too  large,  the  liquid/gas 
interface  disturbance  of  the  jet  can  be  violent  enough  to  invalidate  the  quiescent  interface 
assumption. 

For  the  configuration  of  Figures  1  and  4,  the  significance  of  the  effect  of  the  N2  jet  impinging  on  the 
liquid  surface  can  be  determined  from  an  estimate  of  the  axial  velocity  of  the  orifice  jet  at  the 
liquid/gas  interface.  Such  an  estimate  is  obtained  here  from  the  characteristics  of  an  incompressible 


submerged  jet  (i.e„  analogous  to  the  N2  orifice  jet)  in  an  unconfined  space  (analogous  to  the  gas 
volume  of  the  discharge  vessel).  The  estimate  is  also  useful  for  other  orifice  orientations. 

Let  the  velocity  on  the  jet  axis  a  distance  Zj^j  =  ZDV  -  Z  from  the  oriSce  (i.e.,  at  the  elevation  of 
the  liquid/gas  interface)  be  denoted  by  U  ^  and  assume  a  uniform  velocity,  UQ,  across  the  orifice 
opening 


U0  -  (dMDVN>/dt)/[A0(MHx>N2'/^HT)] 


(13) 


From  [6] 


UjetAJq  =  0.96/[0.066(2Zjet/Do)  +0.29]  for  2ZJET/D0  >  10.  (14) 

When  the  time-dependent  solution  for  MDyN2  ^  available,  Uj^j  can  now  be  estimated  from  Eqs. 
(13)  and  (14)  and  the  significance  of  the  jetting  phenomenon  can  be  assessed. 


SOLVING  THE  MODEL  EQUATIONS 

A  computer  program  was  developed  to  solve  the  initial  value  problem  of  Eqs.  (11)  and  (12) 
corresponding  to  an  arbitrary  choice  of  geometric  parameters,  material  properties,  and  initial 
conditions,  and  to  determine  the  value  of  Uj^  from  Eqs.  (13)  and  (14).  The  method  of  solution  is 
based  on  the  differential  equation  solver  RKQC  presented  in  [7].  The  program  was  used  in  the 
example  calculations  to  follow. 


EXAMPLE  CALCULATIONS 

An  Experimental  Procedure  to  Simulate  the  Discharge  of  Field-Deployed  Vessels 

As  mentioned  in  the  introductory  comments,  the  example  calculations  to  be  presented  here  focus  on 
the  establishment  of  an  experimental  design  and  procedure,  where  vessel  discharge  closely  simulates 
the  discharge  of  a  field-deployed  vessel  while  allowing  for  acquisition  of  data,  including  high  speed 
photography,  to  characterize  the  discharge  process. 

Field-deployed  systems  use  an  explosive  device  to  remove,  "on  demand,"  a  cap  covering  the  exit 
nozzle/orifice  of  a  pre-press urized  discharge  vessel.  Thus,  the  deployed  system  is  characterized  by 
a  Figure  1-type  configuration  with  no  holding  tank  component,  Le.,  with  A 0  -  0. 

The  nature  of  the  experimental  program  at  NIST  precludes  the  use  of  an  explosive  cap  device.  The 
experimental  procedure  to  be  evaluated  involves  a  vessel  discharge  process  which  is  initiated,  instead, 
by  rupture  of  a  nozzle/orifice  diaphragm  cap  at  a  high  cross-diaphragm  pressure  difference.  The 
pressurization  of  the  discharge  vessel  from  PDV1  to  *be  diaphragm  burst  pressure,  Pburst>  “ 
achieved  with  the  use  of  the  holding-tank/orifice-flow  feature  of  Figure  1. 
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As  discussed  earlier,  an  experimental  run  involves  sequential  processes  of  vessel  pressurization  and 
discharge.  Initially,  the  diaphragm  cap  prevents  flow  from  the  vessel  to  the  outside  environment.  At 
t  =  0  the  orifice  connecting  the  vessel  and  the  holding  tank  is  opened  and  pressurization  of  the 
vessel  is  initiated.  At  the  instant  the  pressure  in  the  vessel  reaches  the  diaphragm  burst  pressure,  the 
diaphragm  ruptures.  Nozzle/orifice  flow  and  discharge  of  the  vessel  are  then  initiated.  As  mentioned 
earlier,  the  pressurization  and  discharge  processes  of  interest  here  are  expected  to  occur  over  time 
intervals  of  the  order  of  10s  and  10'2s,  respectively 

Prior  to  any  particular  test  run,  Pburst ’s  onty  known  to  an  accuracy  of  approximately  ten  percent. 
Also,  the  high  speed  camera  to  be  used  to  photograph  the  discharge  process  during  the  test 
procedure  can  record  a  total  time  interval  no  larger  than  the  order  of  a  few  seconds,  including  an 
initial  interval  of  approximately  Is  required  to  bring  the  camera  up  to  its  operating  speed. 


Criteria  for  Experimental  Discharges  to  Closely  Simulate  Field-Deployed  Discharges 

One  would  hope  to  de'ermine  that  an  experimental  discharge  would  closely  simulate  a  field-deployed 
system  discharge  by  obtaining  good  agreement  between  results  of  model  simulations  of  the  two 
processes.  The  basic  criterion  would  be  good  reproduction  of  the  predicted  time-dependent  values 
of  PDV  and 

Since  field-deployed  discharges  involve  no  addition  of  N2  and  no  significant  heat  or  mass  transfer  at 
the  liquid/gas  interface,  an  acceptable  Figure  1-type  experimental  procedure  should  similarly  involve 
discharges  with  no  significant  addition  of  N2.  In  view  of  this,  in  addition  to  the  above  criterion  two 
additional  criteria  must  be  satisfied  if  the  experimental  discharge  is  to  simulate  closely  a  field-deployed 
discharge:  1)  the  experiment  must  include  a  discharge  process  during  which  the  total  mass  of  N2 
delivered  from  the  holding  vessel  to  the  discharge  vessel  is  small  compared  to  the  mass  of  N->  in  the 
discharge  vessel  immediately  prior  to  the  onset  of  discharge,  i.e.. 


require  (MN,  -  MN23)/MN^B  <<  1  for  t  >  tBURST 


(15) 


where  MN,  B  is  the  mass  of  N2  in  the  discharge  vessel  at  tBURST;  and  2)  the  addition  of  N,  during 
the  entire  pressurization/discharge  sequence  does  not  lead  to  gas  flows  in  the  discharge  vessel  gas 
(where  velocities  will  be  of  the  order  of  Ujg-p)  which  are  so  vigorous  as  to  lead  to  significant  heat  or 
mass  transfer  at  the  liquid/gas  interface.  It  is  reasonable  to  expect  that  interface  surface  interactions 
will  not  be  significant  if  UJET  -  dZ/dt  never  exceeds  the  order  of  a  few  m/s.  Thus 


require  Ujgj  -  dZ/dt  <  5  to  lOm/s  for  t  >  0 


Model  Input  Parameters 

Sets  of  parameters  representative  of  a  possible  Figure-l-type  experimental  design  and  associated 
design  procedures  are  presented  in  Table  1.  The  parameters  are  selected  to  simulate  discharge  of 
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a  field-deployed  system  involving  a  0.5xl0'3m3  discharge  vessel  half-filled  with  liquid  CHC1F2  and 
pressurized  with  N2  to  41.37xlOra  (600psi).  Parameters  of  the  field-deployed  system  are  included 
in  Table  1. 

The  parameters  sets  of  Table  1  define  simulations  of  18  different  experimental  systems  and  the  single 
simulation  of  the  field-deployed  system.  The  parameters  sets  were  used  to  define  and  solve  the 
model’s  initial  value  problem.  Solution  results  provide  time-dependent  histories  of  all  model  variables 
and  the  information  necessary  to  determine  whether  or  not  the  criteria  of  Eqs.  (15)  and  (16)  are 
satisfied  during  the  simulation.  Solution  results  will  be  presented  and  discussed  below. 

As  indicated,  the  model  simulations  involve  the  agent  CHCffV  In  solving  Eqs.  (11)  and  (12)  the 
following  thermodynamic  properties  were  used  for  this  and  for  N2. 


Properties  of  Chlorodifluoromethane/CHgF;  (Freon  221  [8], 


Molecular  Weight  =  86.469kg/(kg-mole) 

Density  of  Liquid  =  =  pSAT(294K)  =  1209kg/m3  (17) 

Specific  Heat  (of  gas)  at  Constant  Pressure  (at  300K)  =  CP^G  =  57(103)J/(kg-mole*K) 
PSAr  =  6895(T/R)'C[F  -  (T/R)]E[F  *  (T/HJJ/IFCX/R)]  iqA  -  B/(T/R)  +  D(T/R)pa 

where 


A=29.35754453;  B=3845.193;  C=7.86103122; 
D =0.002190939044;  E=305.826813;  F=686.1 


Properties  of  Nitrogen  (9). 

Molecular  Weight  =  28.013kg/(kg-mole) 

Specific  Heat  at  Constant  Pressure  =  CP^2  (IS) 

=  [7.440 -324(10'3)(T/K)  +  6.400(10-6)(T/K)2  -  1790(10'9)(T/K)3]4186J/(kg-mole  -K) 

where  T  =  Tdv,i  was  used  to  estimate  (a  constant  value  for)  CP(N2. 
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Other  Assumptions  on  the  Model  Parameters 

The  following  completes  the  information  used  to  fully  define  the  model  equations: 

Bow  coefficient  for  the  nozzle/orifice  (flow  of  test  agent  liquid  from  the  discharge  vessel  to  the 
outside  environment!. 


C-D.N  -  0-60 


(19) 


Spinoidal  pressure.  Use  of  Eq.  (8)  and  the  ax  and  o3  components  of  Eqs.  (11)  and  (12)  require 
values  cf  Psp  for  isentropic  paths  from  the  time-dependent  liquid  states  in  the  discharge  vessel.  For 
the  present  calculation  it  is  assumed  that  throughout  the  discharge  process  Psp  <  P^,^  =  1.01x10% 
(14.7psi),  i.e.,  spinoidal-curve  type  instabilities  of  the  metastable  exit  liquid  stream  do  not  occur  during 
the  major  portion  of  the  discharge-process  time  interval  This  assumption  is  consistent  with 
previously  mentioned  recent  Freon  22  discharge  experiments  at  NIST  [5].  However,  as  seen  above 
in  the  discussion  of  C02.  this  assumption  can  not  be  expected  to  hold  in  general  and  actual  estimates 
of  Psp  will  have  to  be  included  in  more  general  applications  of  the  mathematical  model 

Bow  coefficient  for  the  orifice  (flow  of  N-,  from  the  holding  tank  to  the  discharge  vessel!.  CD  0  for 
compressible  flow  through  sharp-edged  orifices  as  a  function  of  cross-orifice  pressure  ratio,  PDv/Pht 
is  provided  in  [10].  Data  points  for  Cp  p(P Qy/P are  plotted  in  Figure  5.  In  the  present  calcula¬ 
tions,  values  of  CD  0  are  estimated  by  linear  interpolation  between  these  data  points. 

Input  value  for  MDVAG.  To  determine  the  constant  mass  of  gaseous  test  agent  in  the  discharge 
vessel  it  is  assumed  that  the  process  of  filling  the  vessel  included  a  time  when  pure  test  agent  filled 
the  entire  vessel  with  gas  and  liquid  volumes  in  thermodynamic  equilibrium  with  each  other.  The 
mass  of  gaseous  agent  at  this  time  is  taken  to  be  In  the  present  calculations  it  is  assumed 

that  at  this  time  of  equilibrium  the  agent  temperature  was  294K.  This  is  the  value  taken  above  for 
the  initial  temperatures  TDV1  and  T^. 

RELATIVELY  LOW-PRESSURE  LARGE-VOLUME  HOLDING  TANK 

This  section  will  present  results  for  the  pressurization/discharge  sequence  using  a  relatively  low- 
pressure  large-volume  holding  tank.  Note  that  experimental  safety  considerations  motivate  use  of 
relatively  low-pressures  in  the  holding  tank,  especially  when  it  is  of  relatively  large  volume. 

Vht  =  2.5xl0'3m3  is  chosen  to  be  ten  times  larger  than  the  initial  gas  volume  in  the  discharge  vessel 
^hxi  =  51.71x10%  (750psi)  is  chosen  to  be  6.89xl05Pa  (IQQpsi)  greater  than  the  marimum  of 
the  diaphragm  burst  pressure  which  can  fall  within  the  range  37.92x10%  <  Pn,m<rr  <  44.82x10% 
(550psi  <  PBURST  <  650psi).  BURST 

The  objective  in  the  present  choice  of  parameters  is  to  provide  for  a  pressurization  process  which  is 
slow  enough  to  minimize  disturbances  of  the  liquid/gas  interface,  Le.  satisfy  Eq.  (16),  but  fast  enough 
to  be  able  to  bracket  clearly  the  time  of  the  subsequent  discharge  process  to  a  known  time  interval 
of  the  order  of  Is.  The  Is  time  interval  is  required  to  guarantee  timely  triggering  of  and 
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photographic  data  acquisition  from  a  high-speed  camera  during  the  discharge  process.  As  mentioned 
above,  the  latter  is  expected  to  occur  over  a  time  interval  of  the  order  of  10'2s. 

The  relatively  large  volume  of  the  holding  tank  would  allow  the  pressurization/discharge  process  to 
be  initiated  either:  1)  from  a  relatively  high-PDV1  state,  somewhat  below  the  minimum  possible 
diaphragm  burst  pressure  and  with  a  relatively  significant  initial  mass  of  N2  in  the  discharge  vessel; 
or  2)  from  a  minimum-PDV1  state  where  there  is  no  N2  in  the  discharge  vessel  (i.e,  the  discharge 
vessel  initially  contains  only  pure  test  agent  at  the  saturation  pressure  PDVj  =  Pv,aT.  where  PsaT  is 
significantly  below  the  Pburst  pressures  under  investigation). 

Solution  res ul's  are  presented  in  Figures  6  -  13  and  these  are  discussed  below. 


Initiating  the  Pressurization/Discharge  Sequence  from  the  Agent’s  Saturation  State;  No  N2  in  the 
Discharge  Vessel  at  t  =  0 

Simulations  were  carried  out  for  experiments  where  the  pressurization/discharge  processes  were 
initiated  at  the  test  agent’s  saturation  state,  i.e„  at  PDV1  =  P$at(Tdv,i  =  294K)  =  9.38xl05Pa 
(136psi),  e.g.,  shortly  after  filling  the  evacuated  vessel  with  the  test  agent.  In  such  experiments  there 
is  no  N2  in  the  discharge  vessel  at  t  =  0,  i.e.,  MDVN->  1  =  0.  Results  of  the  calculations  are  plotted 
in  Figures  6-11. 

The  pressurization  process.  Figure  6  presents  plots  of  PDV  and  during  vessel  pressurization,  0 
<  t  <  tBURST.  Results  are  presented  for  the  three  orifice  diameters,  D0  =  0.005m,  0.001m,  and 
0.0005m. 

As  seen  in  Figure  6,  for  the  smallest  orifice,  D0  =  0.0005m  the  pressurization  process,  when  Pnyi 
-  ?dv  <  PbursT’  takes  4.2s  [for  Pburst  =  37.92xl05Pa  (550psi)]  to  5.8s  [for  Pburst  =  44.82xl(rPa 
(650psi)].  For  the  intermediate-size  orifice,  D0  =  0.001m,  the  pressurization  process  takes  1.0s  (for 
Pburst  =  550psi)  to  1.4s  (for  Pburst  —  650psi).  For  the  largest  orifice,  D0  =  0.005m,  the  entire 
process  takes  place  within  0.1s. 

Ujet  values  corresponding  to  Figure  6  are  presented  in  Figure  7.  For  the  smallest  orifice,  D0  ~ 
0.0005m,  UjgT  starts  out  at  approximately  4m/s  and  drops  to  2  -  3m/s  at  tBURST,  the  time  of 
diaphragm  bursting.  For  D0  =  0.001m,  Uj£T  starts  out  at  approximately  8m/s  and  drops  to 
approximately  4m/s  at  tBURST.  For  D0  =  0.005m,  UJET  starts  out  at  approximately  40m/s  and  drops 
to  as  low  as  approximately  16m/s  [for  Pburst  =  44.82xl0sPa  (650psi)]  at  tBURST. 

Based  on  the  above  results  it  is  concluded  that  during  the  pressurization  process  and  relative  to  the 
criterion  of  Eq.  (16):  1)  use  of  the  large  orifice  would  lead  to  conditions  which  do  not  satisfy  the 
criterion;  2)  use  of  the  intermediate  size  orifice  would  lead  to  results  which  are  barely  acceptable;  and 
3)  use  of  the  small  orifice  would  lead  to  acceptable  performance,  provided  that  pressurization  times 
are  consistent  with  timing  requirements  of  the  photographic  data  acquisition  system. 

The  discharge  process.  Figures  8-11  present  plots  of  Uj^,  dZ/dt,  ZJZV  and  (MN,  -  MNiB)/MN2iB 
during  the  discharge  process.  For  the  limiting  burst  pressures,  Pburst  =  37.92xlCrPa  (550psi)  and 
44.82xl05Pa  (650psi),  results  are  presented  for  the  small-  and  large-orifice  designs. 


Figure  S  presents  results  for  PBURSX  =  3 7.92x10' Pa  (550psi)  and  DQ  =  0.0005m.  Once  the  discharge 
process  is  initiated,  the  liquid  is  seen  to  drop  from  its  initial  elevation  (Z/Zj  =  1)  and  to  be  removed 
entirely  from  the  vessel  (Z/Zj  —  0)  in  approximately  0.024s.  Liquid  discharge  occur;  at  a  relatively 
uniform  rate  with  dZ/dt  dropping  from  approximately  7.0  to  4.5  m/s  during  the  emptying  process. 
During  discharge  UJFT  first  rises  slightly  (because  of  relatively  rapid  reductions  of  PDV,  slow 
reductions  of  Pirr.  and  resulting  increases  in  cross-orifice  pressure  difference)  from  approximately  3 
m  s.  and  then  drops  (because  of  the  significantly  increased  length  along  the  jet  axis  between  the 
source  of  the  jet.  at  the  orifice,  and  the  elevation  at  the  gas/liquid  interface  -  refer  to  Figures  1  and 
2)  to  the  final  value  of  approximately  2m/s.  As  can  be  seen  from  the  plot  of  (MN2  -  MN,  B)/MN7  B. 
during  the  entire  emptying  process  the  total  mass  addition  of  N2  from  the  holding  tank  to  the 
discharge  vessel  is  less  than  one  percent  of  the  initial  N2  mass  in  the  discharge  vessel. 

Figure  9  presents  results  for  PBURST  =  37.92xl05Pa  (550psi)  and  D0  =  0.005m.  Based  on  an 
approximate  extrapolation  from  t  -  tBURST  =  0.020s,  once  the  discharge  process  is  initiated  the  liquid 
is  seen  to  drop  from  its  initial  elevation  and  to  be  removed  entirely  from  the  vessel  in  approximately 
0.023s.  Liquid  discharge  occurs  at  a  relatively  uniform  rate  with  dZ/dt  dropping  from  approximately 
7.0  to  5.5m/s  during  the  emptying  process.  During  discharge  Ujej  first  rises  slightly  from 
approximately  30m's,  and  then  drops  to  the  final  value  of  approximately  20m/s.  As  can  be  seen  from 
the  plot  of  (MN2  -  MN2>B)/Mfj,  B,  during  the  entire  emptying  process  the  total  mass  addition  of  N-, 
from  the  holding  tank  to  the  discharge  vessel  is  significant,  approximately  45  percent  of  the  initial 
mass  in  the  discharge  vessel. 

Figure  10  presents  results  for  FBURST  =  44.82xl05Pa  (650psi)  and  D0  =  0.0005m.  Here  the 
discharge  is  completed  in  0.023s.  Liquid  discharge  occurs  at  a  relatively  uniform  rate  with  dZ/dt 
dropping  from  approximately  7.5  to  5.0m/s  during  the  emptying  process.  During  discharge  UJET  first 
rises  from  approximately  1  to  2m/s,  and  then  drops  back  to  the  final  value  of  approximately  lm/s. 
During  the  discharge  the  total  mass  addition  of  N2  from  the  holding  tank  to  the  discharge  vessel  is 
less  than  one  percent  of  the  initial  N2  mass  in  the  discharge  vessel. 

Figure  1 1  presents  results  for  PBURSX  =  44.82xlOsPa  (650psi)  and  D0  =  0.005m.  From  this  it  is  seen 
that  discharge  occurs  in  approximately  0.021s.  Liquid  discharge  occurs  at  a  relatively  uniform  rate 
with  dZ/dt  dropping  from  approximately  Z5  to  5.5m /s  during  the  emptying  process.  During  discharge 
UJET  first  rises  fro®  approximately  17  to  25m/s,  and  then  drops  to  the  final  value  of  approximately 
20m/s.  During  the  discharge  the  total  mass  addition  of  N2  from  the  holding  tank  to  the  discharge 
vessel  is  approximately  40  percent  of  the  initial  N,  mass  in  the  discharge  vessel 

Based  on  the  above  results  it  is  concluded  that  during  the  discharge  process  and  relative  to  the 
criteria  of  Eqs.  (15)  and  (16):  1)  use  of  the  large  orifice  would  lead  to  conditions  which  violate  both 
of  the  criterion;  2)  use  of  the  small  orifice  would  lead  to  acceptable  performance. 


Initiating  the  Pressurization/Discharge  Sequence  From  a  Slightly  Below  PBURST 

Simulations  were  carried  out  for  experiments  initiated  from  a  PDV1  which  is  slightly  below  PBUrst- 
Since  the  minimum  possible  PBURST  is  37.92xl05Pa  (550psi),  P™,  is  selected  to  be  34.47xl(PPa 
(500psi). 
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The  pressurization  process.  Figure  12  presents  plots  of  PDV  and  during  vessel  pressurization. 
Results  are  presented  for  the  thiee  orifice  diameters,  Dq  =  0.005m,  0.001m,  and  0.0005m. 

As  seen  in  Figure  12.  for  the  smallest  orifice,  D0  =  0.0005m,  and  for  the  diaphragm  burst  pressures 
under  consideration,  the  pressurization  process,  Pov.i  -  Pdv  <  Pburst*  takes  0.6s  [Pburst  = 
37.92xl05Pa  (550psi)]  to  1.9s  [Pburst  =  44.82xl05Pa  (650psi)].  For  the  intermediate-size  orifice,  D0 
=  0.001m,  the  pressurization  process  takes  0.1s  (Pburst  ~  550psi)  to  0.5s  (Pburst  =  650psi).  For 
the  largest  orifice,  D0  =  0.005m.  the  entire  pressurization  process  takes  place  in  less  than  0.02s. 

Ujet  v3lues  corresponding  to  Figure  12  are  presented  in  Figure  13.  For  the  smallest  orifice,  D0  = 
0.0005m,  UJET  starts  out  at  approximately  3.5m/s  and  drops  to  2.0  -  2.5m/s  at  the  time  of  diaphragm 
rupture.  For  D0  =  0.001m.  Ujet  starts  out  at  approximately  7m/s  and  drops  to  approximately  4.5m/s 
at  the  time  of  diaphragm  rupture.  For  D0  =  0.005m,  Ujet  starts  out  at  approximately  33m/s  and 
drops  to  as  lew  as  approximately  21m/s  [for  Pburst  =  44.82xl05Pa  (650psi)]  at  the  time  of  diaphragm 
rupture. 

Based  on  the  above  results  it  is  concluded  that  during  the  pressurization  process  and  relative  to  the 
criterion  of  Eq.  (16),  the  present  results  for  PDV1  =  34.47xl05Pa  (500psi)  are  similar  to  the  earlier 
results  for  PDV1  =  PSAT  =  9.38xl0sPa  (136psi). 

The  discharge  process.  In  the  present  case  where  PDV1  =  34.47xl05Pa  (500psi),  results  for  the 
discharge  process  are  very  similar  to  the  earlier  results  where  PDV1  =  PsaT  =  136psi  and  are 
qualitatively  well  represented  by  the  plots  of  Figures  8-11  and  the  earlier  discussion  of  these. 

Comparing  the  Discharges  of  a  Small-Orifice  Test  Configuration  and  a  Field-Deployed  System 

The  model  equations  were  also  used  to  simulate  the  discharge  of  a  field-deployed  vessel  filled  with 
CHC1F2  and  pressurized  with  N2.  At  the  onset  of  discharge  the  temperature  in  the  vessel  was  taken 
to  be  294K  and,  corresponding  to  the  above  calculations,  the  pressure  in  the  vessel  was  taken  to  be 
41.37xl05Pa  (600psi)  [i.e.,  a  nominal  Pburst  =  41.37xl(rPa  (60Gpsi),  where  simulating  tests  would 
involve  experiments  where  37.92xlOsPa  <  Pburst  -  44.82xl(rPa  (550psi  £  Pburst  -  650psi)].  In 
Figures  14  and  15  calculated  values  of  PDV  and  ZJZ1  for  the  deployed  system  simulation  are 
compared  with  the  values  of  these  variables  obtained  in  the  above  test  simulations  with  D0  = 
0.0005m.  (The  PDV1  =  O^SxlO^a  (136psi)  results  are  presented  in  the  figures  since  these  compare 
somewhat  less  favorably  with  the  field-deployed  system  results  than  do  the  PDV1  =  34.47xl05Pa 
(500psi)  results.] 

As  can  be  seen  in  the  figures,  simulated  field-deployed  system  results  and  simulated  test  results 
compare  favorably. 

Summary  of  Results  of  Simulations  Involving  the  Relatively  Low-Pressure  Large-Volume  Holding 
Tank 

The  above  described  model  simulations  lead  to  the  following  summary  result: 


A  Figure  1-typc  test  configuration  with  a  0.25xl0'3m3  holding  tank  and  a  0.0005m,  and 
possibly  a  O.OOlm-diamctcr  orifice  will  provide  experimental  discharges  which  can  be  expected 
to  simulate  accurately  the  discharge  of  field-deployed  systems.  Of  the  two  orifice  sizes  the 
smaller  orifice  is  preferable.  Use  of  a  0.005m  diameter  orifice  can  not  be  expected  to 
adequately  simulate  field-deployed  system  discharges  under  the  conditions  studied. 

It  is  expected  that  this  result  can  be  extended  to  a  wide  range  of  test  parameters  and  test  agents. 


RELATIVELY  HIGH-PRESSURE  SMALL-VOLUME  HOLDING  TANK 

This  section  will  present  results  for  the  pressurization/discharge  sequence  using  a  relatively  high- 
pressure  small-volume  holding  tank.  VfIj  =  2.5xl0'5m3  is  chosen  to  be  one  tenth  of  the  initial  gas 
volume  in  the  discharge  vessel  and  =  155. 13x10^3  (2250psi)  is  chosen  to  be  large  enough  to 
insure  that  the  amount  of  N2  stored  in  the  holding  tank  will  rupture  the  diaphragm  even  when 
Pburst  is  **s  maximum  value.  44.82xl03Pa  (650psi). 

The  objective  of  choosing  this  combination  of  parameters  is  to  attain  a  pressurization/discharge 
sequence  which  is  so  rapid,  well  within  the  above-mentioned  Is  time  interval,  as  to  guarantee,  without 
any  significant  timing  considerations,  successful  acquisition  of  high-speed  photographic  data  of  the 
discharge  process.  The  feature  of  a  relatively  small  initial  volume  and  mass  of  N2  in  the  holding  tank 
is  only  consistent  with  test  sequences  with  relatively  high  initial  PDV1’s,  values  close  to  the  minimum 
possible  diaphragm  burst  pressure.  Thus,  the  present  test  parameters  would  not  be  expected  to  be 
useful  in  initiating  the  pressurization/discharge  process  from  an  initial  state  where  the  vessel  contained 
only  pure  test  agent  at  its  saturation  pressure.  Here,  PDV1  is  again  chosen  to  be  34.47x10^3 
(500psi),  slightly  below  37.92xl05Pa  (550psi),  the  minimum  value  for  Pburst- 

Solution  results  are  presented  in  Figures  16-24. 


Pressurization  and  Discharge  for  the  Small  Orifice,  D0  =  0.0005m 

Simulation  results  for  pressurization/discharge  test  using  the  small  orifice  design  are  presented  in 
Figures  16  -  20. 

F°r  ^burst  =  37.92xl05Pa  (550psi),  plots  of  PDy  and  P^-  are  presented  in  Figures  16  and  plots  of 
UjeT’  dZ/db  Z/Z,,  and  (MN2  -  MN2iB)/MN2>B  in  Figure  17.  Note  that  the  latter  variable  is  meaningful 
only  during  discharge  when  t  >  tBURST.  As  can  be  seen  in  these  figures,  for  this  case  = 

0.173s,  i.e.,  this  is  the  time  it  takes  for  PDV  to  rise  from  PDV1  =  34.47xl05Pa  (500psi)  to  the 
diapnragm  rupture  pressure.  From  Figure  17  it  is  seen  that  the  l’iquid  is  nearly  all  discharged  from 
the  vessel  (i.e.,  Z/Zj  is  approaching  zero)  at  approximately  t  =  0.2s.  From  Figure  17  it  is  also  seen 
that  during  discharge:  the  rate  of  liquid  outflow  is  reduced  to  about  half  of  its  original  value,  with 
dZ/dt  being  reduced  from  7.0m/s  to  4.5m/s;  is  reduced  from  an  initial  pre-discharge  value  of 
approximately  4m/s  to  approximately  2m/s;  and  (M^2  -  M.n2,b)/Mn2,8  rises  from  0  to  approximately 
0.02,  i.e.  at  the  end  of  the  discharge  the  mass  of  N2  in  the  vessel  is  increased  from  its  value  at  tBtJRST 
by  a  factor  of  approximately  0.02.  The  latter  two  results  indicate  that  the  criteria  of  Eqs.  (15)  and 
(14)  are  both  satisfied  and.  in  this  sense,  as  with  the  relatively  large  holding  tank,  the  small  orifice 
design  can  be  expected  to  adequately  reproduce  field-deployed  discharges. 


For  Pburst  =  44.82xl05Pa  (650psi),  plots  of  PDV  and  Pm  are  presented  in  Figures  18  and  plots  of 
Ujet  and  dZ/dt  are  presented  in  Figure  19.  Plots  of  UJET,  dZ/dt,  ZfZv  and  (MN2  -  MN2ib)/Mn2,B 
for  the  times  of  the  discharge  process  are  highlighted  in  Figure  20.  For  this  case  tBURST  =  0.94s. 
From  Figure  18  it  is  seen  that  the  Pm  has  almost  been  reduced  to  Pny  at  the  time  that  discharge 
is  initiated.  Indeed,  if  the  P^ri  had  been  much  less  than  155.13xlO*Pa  (2250psi),  PDV  would  not 
have  risen  to  PguRST  and  the  discharge  process  would  never  have  occurred!  From  Figure  19  it  is 
seen  that  during  the  pressurization  process  UJET  is  reduced  from  approximately  4.0m/s  to  1.5m/s.  The 
liquid  is  nearly  all  discharged  from  the  vessel  at  0.962s.  During  discharge:  dZ/dt  is  reduced  from 
7.5m/s  to  5.0m/s;  UJET  rises  to  approximately  2.5m/s  from  its  initial  value  of  1.5m/s  and  finally  drops 
to  approximately  2m/s;  and  (MN2  -  MN2  B)/MN2  B  rises  from  0  to  approximately  0.005.  The  criteria 
of  Fqs  (15)  and  (16)  are  again  satisfied. 

As  noted,  for  both  PBURSt  =  37.92xl05Pa  (550psi)  and  44.82xl0sPa  (650psi)  the  above  results 
indicate  that  both  criteria  of  Eq.  (15)  and  (16)  are  satisfied  with  the  small-holding-tank/small-orifice 
design  and  that  in  this  sense  this  design  can  be  expected  to  adequately  reproduce  field-deployed 
discharges.  Flowever,  for  the  following  three  reasons  the  large-holding-tank  design  is  significantly 
more  robust  than  the  small-holding-tank  design: 

1.  For  a  given  maximum  value  for  PBURST,  and  discharge  vessel  volume,  the  minimum 
acceptable  value  for  Pht,i  (approximately  lSS.DxlO^a  (2250psi)  for  the  present  test 
parameters]  is  very  sensitive  to  the  original  amount  of  liquid  in  the  vessel.  For  example,  in 
the  present  simulations  the  vessel  was  assumed  to  be  one-half  full;  but  if  it  were  only  one- 
quarter  full,  the  minimum  Pprx;i  that  would  lead  to  rupture  of  the  44.82xl05Pa  (650psi) 
diaphragm  would  have  been  approximately  206.84x1 (^Pa  (3000psi). 

2.  The  small-volume-holding  tank  configuration  does  not  allow  for  the  option  of  initiating  tests 
from  relatively-low  values  of  PDVj,  e.g.  immediately  after  filling  the  vessel  with  the  test  agent, 
when  PDV1  =  Psat-  To  d°  so  would  once  again  require  highly  variable  and  unreasonably 
large  values  of  Phxi- 

3.  The  large-holding-tank  design  requires  P^m  values  only  slightly  greater  than  the  maximum 
PfiURST  value  ancl  it  suffers  from  neither  of  the  above  shortcomings. 


Pressurization  and  Discharge  for  the  Large  Orifice,  D0  =  0.005m 

Simulation  results  for  pressurization/discharge  test  using  the  large-orifice  design  are  presented  in 
Figures  21  -  24. 

For  Pburst  =  37.92xl05Pa  (550psi),  plots  of  PDV  and  P^  are  presented  in  Figures  21  and  plots  of 
UJET,  dZ/dt,  Z/Zj,  and  (MN2  -  MN2^B)/MN2jB  in  Figure  22.  Note  that  the  latter  variable  is  meaningful 
only  during  discharge  when  t  >  tBURST.  As  can  be  seen  in  these  figures,  for  this  case  tBURST  = 
0.00173s,  i.e.,  this  is  the  time  it  takes  for  PDV  to  rise  from  PDV1  =  34.47xl05Pa  (500psi)  to  the 
diaphragm  rupture  pressure.  From  Figure  21  is  seen  that  the  P^  has  almost  been  reduced  to  PDV 
at  t  =  0.015s,  at  which  time  approximately  40  percent  of  the  original  amount  of  liquid  is  still  in  the 
vessel,  i.e.  Z/Zj  =  0.4  (Figure  22).  Also,  from  Figure  22  it  is  seen  that  the  liquid  is  nearly  all 
discharged  from  the  vessel  (i.e.,  ZJZ1  is  approaching  zero)  at  approximately  t  =  0.023s.  From  Figure 
22  it  is  seen  that  during  discharge:  the  rate  of  liquid  outflow  is  reduced  to  about  half  of  its  original 
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value,  with  dZ/dt  going  from  8m/s  to  4m/s;  UJET  is  always  dropping,  being  reduced  from  an  initially 
large  pre -discharge  value  of  approximately  40m/s  to  approximately  4m/s;  and  (MN2  -  MN2  B)/MN2  B 
rises  from  0  to  0.35,  i.e.  at  the  end  of  the  discharge  the  mass  of  N2  in  the  vessel  is  increased  from  its 
value  at  tBURST  by  a  factor  of  approximately  0.35.  The  latter  two  results  indicate  that  the  criteria  of 
Eqs.  (15)  and  (16)  are  not  satisfied  and,  as  when  used  in  conjunction  with  a  relatively  large  holding 
tank,  the  large  orifice  design  is  not  expected  to  adequately  reproduce  field-deployed  discharges. 

For  PBurst  =  44.82xl05Pa  (650psi),  plots  of  PDV  and  Pm  are  presented  in  Figures  23  and  plots  of 
UJET,  dZ/dt,  Z/Zj,  and  (MN2  -  MN2B)/MN,  B  in  Figure  24.  The  results  are  qualitatively  similar  to 
those  discussed  in  the  last  paragraph  for’  PBURST  =  37.92xl05Pa  (550psi).  The  quantitative 
differences  can  be  easily  determined  from  the  figures.  For  this  case  tBURST  =  0.0094s.  From  Figure 
23  is  seen  that  the  P^-  has  almost  been  reduced  to  PDV  at  the  time  that  discharge  is  initiated. 
Indeed,  if  the  P^xi  had  been  much  less  than  155.13xl0ra  (2250psi),  PDV  would  not  have  risen  to 
Pburst  and  discharge  process  would  never  have  occurred!  The  liquid  is  nearly  all  discharged 
from  the  vessel  at  0.03s.  During  pressurization  process  UJET  drops  from  approximately  40m/s  to 
15m/s.  During  discharge:  dZ/dt  is  reduced  from  7m/s  to  5m/s;  UJET  is  reduced  from  15m/s  to  3.5m/s; 
and  (MN2  -  MN2  B)/MN,  B  rises  riom  0  to  0.1.  Although  the  criterion  of  Eq.  (15)  is  satisfied,  that  of 
Eq.  (16)  is  not. 


SUMMARY  AND  CONCLUSIONS 

A  mathematical  model  and  associated  computer  program  were  developed  to  simulate  the  discharge 
of  fire  extinguishment  agents  from  N2-pressurized  vessels.  The  model  is  expected  to  have  three 
applications.  First,  to  establish  an  experimental  procedure  which  both  1)  simulates  de-pressurization 
of  a  field-deployed  discharge  vessel  and  2)  allows  for  acquisition  of  data,  including  high  speed 
photography,  to  chracterize  adequately  the  discharge  process.  Second,  to  evaluate  the  discharge 
characteristics  of  a  wide  range  of  altemative-agent/pressure-vessel  configurations,  thereby  extending 
the  slow  and  relatively  costly  experimental  method  of  making  such  evaluations.  Finally,  to  predict 
exit  flow  conditions  to  be  used  to  solve  the  problem  of  agent  dispersal  outside  of  the  discharge  vessel. 

The  model  is  based  on  the  experimental  configuration  depicted  in  Figure  1,  and  the  solution  method 
is  capable  of  treating  arbitrary  choices  of  geometric  parameters,  material  properties,  and  initial 
conditions. 

The  model  was  used  in  example  calculations  which  address  the  first  of  the  applications,  viz.,  the 
establishment  of  a  robust  experimental  design  and  procedure  that  would  1)  simulate  the  discharge 
of  field-deployed  fire  extinguishment  systems  useable  in  US  Air  Force  aircraft  and  2)  meet  additional 
experimental  constraints  consistent  with  a  program  at  NIST  in  support  of  the  US  Air  Force  fire  safety 
issues. 

In  the  example  calculations,  the  fixed  parameters  which  characterized  the  system  were  a  half-liter 
discharge  vessel  with  a  circular  exit  nozzle/orifice  of  diameter  0.01905m  half-filled  with  Freon  22  in 
the  liquid  phase  at  "room  temperature,"  294K. 

The  deployed  system  was  assumed  to  be  pressurized  with  N2  to  a  total  pressure  of  41.37xl05Pa 
(600psi),  where  the  exit  nozzle/orifice  is  closed  with  an  explosive  cap  which  is  released  "on  demand." 


18 


In  the  experimental  system,  the  exit  nozzle/orifice  is  capped  with  a  disk  which  ruptures  at  an  a  priori- 
unknown  pressure,  between  37.92xl05Pa  (550psi)  and  44.82xl05Pa  (650psi). 

The  experimental  system,  is  brought  to  the  rupture/burst  pressure,  Pburst*  fr°m  an  initial  pressure 
PDvi  by  a  process  of  pressurization  by  flow  from  the  holding  tank  which  is  initially  filled  with  N2  at 
some  pressure,  P^j,  higher  than  PBUrST-  The  holding  tank  communicates  with  the  discharge  vessel 
through  a  circular  orifice  of  diameter  D0. 

The  model  was  used  to  simulate  the  discharge  of  the  deployed  system  and  of  the  experimental  system. 
Experimental  parameters  that  were  varied  were  PDV1  J9.38xlOsPa  (136psi)  or  34.47xl05Pa  (500psi)]; 
Pht.i  anc*  lbe  volume  of  the  holding  tank  [51.71xlCrPa  (750psi)  and  2.5xl0'3m3,  or  lSS.lSxlC^Pa 
(2250psi)  and  2.5xl0‘sm3];  and  D0  (0.005m,  0.001m,  or  0.0005m). 

Results  of  the  simulations  were  presented  in  Figures  6  -  24. 

Based  on  the  model  calculations  it  was  determined  that  with  PDV1  =  9.38x1 05Pa  (136psi)  or 
34.47xl05Pa  (500psi)  an  experimental  configuration/design  using  the  relatively  low-pressure 
[SlJlxlO^a  (750psi)],  large-volume  (2.5xl0'3nri),  holding  tank  and  the  smallest-diameter  (0.0005m) 
orifice  would  accurately  simulate  the  discharge  of  the  field-deployed  system.  Furthermore,  it  is 
expected  that  such  a  design  is  robust  in  the  sense  that  it  would  also  simulate  well  the  field-deployed 
system  even  when  extended  to  a  range  of  parameters  and  agent  materials  well  beyond  the  scope  of 
the  present  calculations. 

Model  calculations  also  indicated  that  an  experimental  design  which  uses  the  large-diameter  orifice 
(0.005m)  and/or  the  relatively-high  pressure  small-volume  holding  tank  [lSS.lSxlO^a  (2250psi)  and 
2.5xl0‘5m3]  would  often  not  provide  accurate  simulations  of  field-deployed  discharges  and  would 
never  be  expected  to  be  "robust." 

Validation  of  the  model  requires  experimental  confirmation  of  model  predictions. 
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3  2 


NOMENCLATURE 


adv 

cross-section  area  of  discharge  vessel 

An,  Aq 

area  of  discharge  nozzle/orifice,  orifice  from  holding  tank 

Cd,n>  ^d.o 

flow  coefficient  of  discharge  nozzle/orifice,  orifice  from  holding  tank 

CP,AG’  Cp^J2 

specific  beat  at  constant  pressure  of  agent  gas,  N2 

CyAG>  ^N2 

specific  heat  at  constant  volume  of  agent  gas,  N2 

Do 

diameter  of  orifice  connecting  discharge  vessel  and  holding  tank 

f1?  f2,  f3,  f4 

functions  of  Xj  and  parameters,  Eq.  (A-3) 

mdv,ag>  ^dyal 

mass  of  gaseous,  liquid  agent  in  discharge  vessel 

MDYN2>  MHTN2 

mass  of  N2  in  discharge  vessel,  holding  tank 

MDV,N2,1>  MDYAL1> 
^PTCN2,1 

MDVN2>  MDYAL>  MHTN2  at  1  ~  0 

MN2 

total  mass  of  N2  in  discharge  vessel  and  holding  tank,  see  Eq.  (3) 

MN23 

mass  of  N2  in  discharge  vessel  when  the  diaphram  ruptures 

P 

pressure 

PAMB 

P  of  ambient  environment 

PBURST 

P  in  discharge  vessel  when  the  diaphram  ruptures 

PDV»  PHT 

P  in  discharge  vessl,  holding  tank 

PDV,1>  PHX1 

Pdv>  Pht  at  t  =  0 

PSAT 

P  when  liquid  agent  is  at  saturation 

Psp 

P  on  agent  spinoidal  curve 

Rag.  rN2 

gas  constants  for  gaseous  agent,  N2 

T 

temperature 

Tal 

T  of  liquid  agent 

tdv.  tht 

T  of  gaseous  agent  in  discharge  vessel,  N2  in  holding  tank 

tdv,i>  tht,i 

1 DV’  "^HT  at  1  “  0 

t 

time  from  begining  of  experiment  or  from  onset  of  discharge  of  field-deployed 
system 

lBURST 

t  when  diaphram  ruptures 

lD 

t  at  completion  of  liquid  discharge  or  when  PDV  reaches  PSAj  during 
discharge 

UJET  U0 

axial  velocity  of  orifice  jet  upstream  of  gas/liquid  interface 

VHT 

volume  of  holding  tank 

Xj.  x2,  x3 

variables  of  dimensionless  initial  value  problem,  Eq.  (9) 

Z 

elevation  of  gas/liquid  interface  above  bottom  of  discharge  vessel 

zi 

Z  at  t  =  0 

ZDV 

length  of  discharge  vessel 

ZJET 

ZDV  '  z 

*  vag 

ratio  of  specific  heats  for  N2,  gaseous  agent 

C 

dimensionless  Z,  Eq  (A-l) 

^DV 

dimensionless  ZDV,  Eq.  (A-4) 

^DV> ^HT 

dimensionless  TDV,  Tjjj,  Eq.  (A-l) 

^DV,1>  0HX1 

dimensionless  TDV1,  T^j,  Eq.  (A-l) 

Pal 

density  of  liquid  agent 

PSAT 

density  of  saturated  liquid  agent 

4 

dimensionsless  parameters,  i  =  1  to  5,  Eqs.  (10)  and  (A-4) 

^DYN2>  PHTN2 

dimensionless  MDVN2,  Eq.  (A-l) 

PDV,N2,1>  ^HXN2,1 

dimensionless  MDVjN2j1,  Eq.  (A-l) 

^AMB 

dimensionless  P^^g,  Eq.  (A-4) 

^HT-  ^DV-  ^HXl 

dimensionless  Pj^,  PDV,  Phxi,  Eq-  (A-l) 

dimensionless  PSP.  Eq.  (A-4) 

functions  of  Xj  and  parameters,  i  =  1,  2,  3,  Eqs.  (11)  and  (12) 
dimensionless  t,  Eqs.  (9)  and  (A-l) 
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APPENDIX:  THE  DIMENSIONLESS  INITIAL  VALUE  PROBLEM 
Define  the  following  dimensionless  variables  and  their  initial  values: 

T  2  t<-'D,NjdLN[2PDV,l/(^DV^L,lADvZ1)]1/2 

77 HT  -  PhT^Dyri  TrHT^T  =  °)  -  ^HTl  =  PhTi/PdV.I  >  1 

■^dv  =  Pdv^P bv,i»  7rDv(T  =  0)  =  l 

PHXN2  2  MHi,N2'7viN2;  PhWt  =  °)  =  MHXN2,1/MN2  -  ^HTN2,1  *  1  *  PDV.N2.1  5  1 
PDV.N2  =  ^DV,N2'^N2>  PDV,N2(t  =  0)  =  ^DV,N2,1^N2  2  ^DyN2,l  =  1  *  11HTN2,1  5  1 
PDyAL  2  MDyAL^DyALl  =  MDyAL/(>ALADVZl);  PdV,Al(t  *  °)  =  1  (A-l) 

®ht  =  "Pht^N2^N2^0^  htPdv.i)  =  CPht^hi;iXPhti^Pdv,i)/,(^htn2,i^n2); 

^Ht(t  =  0)  s  ®HX1  =  (PhT1^PdV,.’^(^HTN2,1^N2)  >  1 
®DV  =  TdvMdV^G^AgAZ1ADvPdV,i)  =  ^DV^DyAG^AC^AL^^DyALlPoyi)’ 

®Dv(t  =  0)  =  ^DV.l  ~  TDV,iMDV>agRAg/(ZiAdvP DV1) 

$  =  Z /Zl5  £(r  =  0)  =  1 

Define  the  new  variable  xx  and  introduce  x2  and  x3  as  new  designations  for  |inv.N2  anc*  ^dval* 
respectively: 

Xl(T)  2  K1  +  ^DVJ^C1  +  ^1  lIDyN2)]7rDv(^DV  *  MdYAl)  +  *2*40  *  PDyN2)ri 

Xl(°)  =  10  +  *3PDyN2,lV0  +  *lPDV,N2,l)K?DV  -  !)  +  W1  -  PDYN2,l)y 

(A-2) 

x2(t)  “  PDMN2’  xz(Q)  —  PDyN2,l 

x3(t)  2  PdV,AL>  *3(0)  =  PDyAL,! 

Define  also  the  following  functions  of  the  Xj: 

^1  =  txl  *  *2*40  *  ^O  +  *1X2M0  +  *3x2)0dV  *  *s)]; 

,  „  (A-3) 

f2  =  A2(l  -  Xj)  ;  fj  =  A^l  -  x^7;  f4  =  A2(l  -  x^/fjj 
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Asf3{y[2/(y  +  l)](r  +  W  •  D/f2}^; 


iff4  2s  [(y  +  l)/?]^'1)! 


h(*i)  = 


«*.)■ 


where 


A5f3(l/f4)^{2y[l  -  (l/f4)^  *  X^/I(y  -  l)y}1/2;  if f4  <  [(y  +  W\Wr '1)] 


(fi  -  ^amb) 


1/2. 


if  7 r 


AMB 


17, 


SP 


(Ar3  oont’d) 


(fl  *  '"’SP)1'2:  *  ’’’AMB  <  ^SP 


*1  =  (MN2/MDViAG)(RN2/RAG);  h.  -  7rHXl^^HXN2,l)y* 

A3  =  (CV^2/C^G)(MN2/MDV^G^  =  ^l^AG  *  IVO*  *  1)> 
a4  =  [(yAG  -  mr  -  iJItVHT/CZiADv)]; 

A5  =  (CD>c>CDiN)(A<yAN)[(MDV>AL/MN2)(ADVZ1/Vm.)^]1/2; 

?Dv  =  Zdv^i  ■  constant  >  1;  yAG  =  C^q/C^^q  ~  rag/°v^G  +  V> 
■"■AMB  =  ^AMB^DV,!’  ^SP  =  PsP^Dyi 

Then  the  Xj  can  be  determined  from  the  solution  of  the  following  initial  value  problem 


dx3/dr  =  -  (yAG  -  l)f2^;  dxj/dr  =  f5;  dxj/dr  =  -  ^ 

xl(0)  =  l(f  +  ^^DyNilVC^  +  ^lHDV,N2,l)K£DV  '  f)  +  ^4(1  *  l1DV>N2,l)y ’  (A*5) 

x2(°)  =  l‘DV^J2a;  x3(°)  =  ►‘DYAL.l 

and  the  solutions  for  the  original  dimensionless  variables  of  Eqs.  (A-l)  can  be  retrieved  from 


1»DV  ~  fl?  ^HT  ~  ^3>  1*HXN2  ~  f  *  x2>  ^DYN2  =  x2>  llDVrAL  “  £  ~  x3> 
<?HT  =  f2;  ^DV  =  (£dV  *  x3)V(f  +  ^lx2) 


(A-6) 


From  the  latter  solutions  the  dimensional  variables  can  be  determined  finally  from  the  definitions  of 
Eqs(A-l). 

The  problem  formulation  of  Eqs.  (A-5)  and  its  solution  are  valid  up  to  the  value  of  r  where  =  0 
(the  last  of  the  liquid  is  driven  from  the  discharge  vessel)  or  f,  is  reduced  to  Psat'^dvi  (the  liquid 
in  the  discharge  vessel  is  subjected  to  a  condition  where  an  equilibrium  state  would  lead  to  flashing). 
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agent 

material 

Chlorodinuoromcthane/CHClFj  (Freon  22) 

TDV,1  ’ 

THT1  “ 

294K 

diameter  of 
discharge 
vessel  * 

DDV  = 

0.05m 

volume  of 
discharge 
vessel  = 

VDV  = 

0.5xl0'3m3 

initial 
volume  of 
liquid 
agent  in 
discharge 
vessel  * 

Vov/2  =  0.25xlG'3mJ 

DN  = 

0.01905m 

VHT  = 

2^xl0'3m3 
(large  holding  tank) 

rsxio^ma 
(small  holding  tank) 

field  simu¬ 
lation 
(no  hold¬ 
ing  tank) 

PHT.l  = 

51.71x1 0*pa 
(750psi) 

155.13xl05pa 
(2250 psi) 

* 

PDV,i  = 

34.47xlOspa 

(SOOpsi) 

(Figs.  12  and  13) 

9.38xl05pa 
(136  psi) 

a  PSAT <Tuvyl) 

(Figs.  6  and  7) 

34.47xi0*pa 

(SOOpsi) 

41.37* 

105pa 

(600psi) 

PBURST  = 

37.92x 

105pa 

(550psi) 

44.82* 

105pa 

(650psi) 

37.92x 

lO^pa 

(550psi) 

44.82x 
105pa 
(650 psi) 

37.92x 

lO5?* 

(550psi) 

44£2x 

105pa 

(650psi) 

41.37* 
105pa 
(600psi) 
(Figs.  14 
and  15) 

»0* 

0.0005m, 
0.001m, 
and  0.005m 

0.0005m, 
0.001m, 
and  0.005m 

0.0005m 

(Figs-S, 

14  and  15), 
0.001m, 
and  0.005m 
(Fig.  9) 

0.0005m 
(Figs.  10, 
14,  and 

15), 

0.001m, 
and  0.005m 
(Fig.ll) 

0.0005m 
(Fig  16  and 
17)* 

0.001m, 
and  0.005m 
i  (Figs.  21 
and  22) 

0.0005m 
(Figs.  IS, 

19,  and 

20), 

0.001m, 
and 
0.005m 
(Figs.  23 
and  24) 

:'v-I 


I 


r  v 


*■:■■■: 


Mm 

I . 

Iv\-  • 

P-1 


®pi 


Table  1.  Characteristics  of  simulated  discharges  -  model  input  parameters 


P-V  diagram  for  a  generic  material;  the  spinoidal  curve  and  regions  of  metastable 
liquid  and  vapor. 
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CRITICAL 
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curve 
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Plots  of  Utj-j- for  OS  t  £  tBURST  and  D0  =  0.005m,  0.001m,  and  0.0005m;  where  V 
=  2.5xl0'3m  ,  TDV1  =  Tot.  =  294K,  PDV1  =  Pc^  =  9J8xlOsPa  (136psi),  PH 
=  51.7Ixl0sPa  (750psi),  P^  £  44.82xYo5Pa  (650psi).  ” 


0  0.005  0.01  0.015  0.02  0.025 


*  *  'burst  M 


Plots  of  Ujg-p  dZ/dt,  ZJZV  and  (MN2  -  MN2?b)/^N23  ^  *burst>  where  Vjjj  = 

2.5x10  3m  ,  Tjjyj  =  Tjn^  =  294K,  Pjy^  =  ^sat  ~  9.38xl05Pa  (136psi),  P ht;*  = 
51.71xlOsPa  (750psi),  PBURST  =  37.92xlOsPa  (550psi)  and  D0  =  0.0005m. 


0  0.005  0.01  0.015  0.02  0.025 


*  "  *BURST 


Plots  of  Uj£j,  dZ/dt,  Z/Zj,  and  (MN2  *  MN2jB)/MN23  for  t  2  tBURST;  where  V^j.  — 
2.5x1 0'3m,  TDyi  =  THX1  =  294K,  PDyi  =  PSAX  =  WSxlCPPa  (13dpsi),  P^  = 
51.71xlOsPa  (750psi),  PBURSX  =  37.92xlOsPa  (550psi)  and  D0  =  0.005m. 


UJET  [m/s],  dZ/dt  [m/s],  Z/Z,,  (M 


Figure  10.  Plots  of  Ujj-p  dZ/dt,  Z/Z^  and  (M^  *  Mn2,bV^N23  for  t  ^  ^burstp  ^ht  " 
2.5xl0'3m  ,  Tdv>1  =  Tmi  =  294K,  Pm;i  =  PSAT  =  ^Sxl^Pa  (136psi),  Pm;i  = 
SlJlxl^Pa  (750psi),  PBURSX  =  44.82xlOsPa  (650psi)  and  D0  =  0.0005m. 


/  (MN2  •  Mn2,b)/Mn2,B 


0.005  0.01  0.015  0.02  0.025 


BURST 


Figure  11. 


Plots  of  Un-p  dZ/dt,  2VZj,  and  (MN2  *  B  for  t  2  tBURST;  where  Vjjj  = 

2-5X10-W,  Tmi  =  Tjn-!  =  294 K,  Pmi  =  PM  =  PJSxl^Pa  (136psi),  P™  = 
->1.71xl(rPa  (750psi),  P rijrct  "  44.82xlO^Pa  (650psi)  and  —  0.005m. 
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Figure  12.  Plots  of  PDV  and  Pjjj  for  0<  t  S  tBURST  and  D0  =  0.005m,  0.001m,  and  0.0005m; 

where  Vm  =  2.5xl04m3,  =  T^-j  =  294K,  P^  =  34.47xlOsPa  (500psi),  P^, 

=  51.71xlOsPa  (750psi),  P™^-  £  44.82xlOsPa  (650psi). 


0.005m 


0.001m 

Do=0.0005m 


10'1 


t<t 


BURST 


Figure  13.  Plots  of  U 


Plots  of  Uj£j  for  0s  t  S  and  Dq  —  0.005m,  0.001m,  and  0.0005m.;  where 

VHT  =  2.5x1 0‘3m3,  Tmi  m  Tjnj!  =  294K,  PDV1  =  34.47xl05Pa  (500psi),  Pm-,  = 
51.71xl(rPB  (750psi),  PBURSX  £  44.82xlOsPa  (650psi). 
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Figure  14.  Plots  of  PDV  during  discharge  of:  1)  field-deployed  vessel  at  Pburst  =  41.37xlOsPa 
(600psi)  and  test  configuration  vessel  [V^  =  2.5xl0*3m3,  D0  =  0.0005m,  * 

psat(294K)  =  9-38xlO?Pa  (136psi)]  at  2)  PBURCT  =  37.92xlOsPa  (550psi)  and  3) 
44.82xlOsPa  (650psi). 


0  0.005  0.01  0.015  0.02  0.025 


* '  'burst  Is] 


Figure  15.  Plots  of  Z/Zt  during  discharge  of:  1)  field-deployed  vessel  at  Pburst  ~  41-37xlOsPa 
(600psi)  and  test  configuration  vessel  [Vj^j.  =  2.510‘3m3,  D0  =  0.0005m,  PDV,I  = 
Psat(294k>  =  9J8xl(PPa  (136psi)]  at  2)  PBURCT  =  37.92x1 O^Pa  (550osi)  and  3) 
44.82xlOsPa  (650psi). 
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Figure  16.  Plots  of  Pqv  and  ] 
=  Thci  =  294K, 
^BURST  ~  37.92x1 


t[s] 


tBURST=0.173S 


Plots  of  UJET,  dZ/dt,  and  Z/Zl,  for  t  &  0  and  dZ/dt,  Z/Zj,  and  of  (MN2  -  MN2jB)/ 
M.N2.B  for  *  *  ‘burst?  rherc  VHT  =  2.5x1 0  sm3,  D0  =  0.0005m,  T-,,  *  Tor.  = 
294K,  Pjjv!  =  34.47xlOsPa  (SOOpsi),  Pin;i  =  lSS.UxlO5!*  (2250psl),’  and  PBm,cr 
=  S^Zx!©5!*  (550psi).  BWST 


dZ/Dt 


t  [s]  tBURST=0*940$ 


Figure  19.  Plots  of  UJET  and  dZ/dt  for  t  &  0  and  D0  =  0.0005m;  where  =  2.5xl0‘sm3, 
TDMl  ~  THXl  =  294K»  PDVL  =  34.47xl05Pa  (500psi),  =  155.13X10sPa 

(2250psi),  PBURST  =  44.82xl<rPa  (650psi). 
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0.945  0.95  0.955  0.96 


tBURSTss0.940$ 


Figure  20. 


Plots  of  Ujet>  dZ/dt,  Z/Zj,  and  (M^  *  ^n2,b)^Mn2,b  ^or  *  ^  ^burst  an^  ~ 
0.0005m;  where  VHT  =  2.5xlOsm3,  TDV1  =  Tin;i  =  294K,  PDV,  =  34.47xlOsPa 
(SOOpsi),  Pm-,  =  155.13xlOsPa  (2250psi),  and  PBURST  =  44.82xlOsPa  (650psi). 


tBURST=0.00173S  t  [s] 


Figure  21.  Plots  of  Ppy  and  Pirr  for  t  a  0  and  D0  =  0.005m;  where  V^.  =  2.5x1 0'5m3, 

=  THT1  =  294K,  Pny!  =  34.47xlOsPa  (SOOpsl),  =  lSS.LJxK^Pa  (2250psi), 
P BURST  =  37.92xl05Pa  (550psi). 


aZNVM/(BZNW  -  2NIAI)  ‘‘z/Z  ‘[s/uj]  ip/zp  ‘[s/ui]  13rn 


Figure  22.  Plots  of  UJET>  dZ/dt,  and  Z/Zl,  for  t  2j  0  and  dZ/dt,  ZJZV  and  of  (MN2  -  M N2jB)/ 
MN2>b  f°r  *  ^  *BURSl5"he"  ^HT  ~  2.5x10  sm3,  Dq  =  0.005m,  Tpyj  —  Tjjjjj  = 
294K,  PDV1  =  34.47xlOsPa  (SOOpsi),  P^  =  155.13xlOsPa  (2250psi),  and  PBURST 
=  37.92xlOsPa  (550psi). 


Figure  23.  Plots  of  PDV  and  for  t  a  0  and  D0  =  0.005m;  where  VHT  =  2.5xl0*sm3,  TDV1 
=  Vi  =  294K,  Ppyj  =  34.47xlOsPa  (SOOpsi),  =  155.13xlOsPa  (2250psi), 
PBURST  *  44.82xlOsPa  (650psi). 
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Figure  24.  Plots  of  Uj^p  dZ/dt,  and  Z/Zl,  for  t  &  0  and  dZtft,  ZJZV  and  of  (MN2  -  MN2jB)/ 
MN2,b  for  *  *burst5  w^ere  ^ht  =  2.5xlO  sm3,  D0  =  0.005m,  TDyX  =  Tj^j  = 
294K,  P^  =  34.47xlOsPa  (500psi),  P^  =  155.13xlOsPa  (2250psi),  and  PBURST 
=  44.82xlOsPa  (650psi). 


